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Abstract 

Following recent work of Chernov, Markarian, and Zhang, it is known that 
the billiard map for dispersing billiards with zero angle cusps has slow decay 
of correlations with rate 1/n. Since the collisions inside a cusp occur in quick 
succession, it is reasonable to expect a much faster decay rate in continuous 
time. In this paper we prove that the flow is rapid mixing: correlations decay 
faster than any polynomial rate. A consequence is that the flow admits strong 
statistical properties such as the almost sure invariance principle, even though 
the billiard map does not. 

The techniques in this paper yield new results for other standard examples 
in planar billiards, including Bunimovich flowers and stadia. 

1 Introduction 

Lorentz gas models (and the associated discrete time billiard maps) are an important 
class of examples in mathematical physics. Their systematic study from the viewpoint 
of smooth ergodic theory was begun by Sinai [30] in the 1970's. Sinai focused on planar 
periodic Lorentz gases with convex obstacles (dispersing billiards) proving properties 
such as uniform hyperbolicity and ergodicity. The situation is analogous to the much- 
studied case of Anosov and Axiom A systems for which advanced statistical properties 
such as decay of correlations and central limit theorems are by now classical [2l [28| [3Tj . 
In particular, if A C M is a mixing hyperbolic basic set with Gibbs measure \x for an 
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Axiom A diffeomorphism / : M — > M, then for all Holder observables v , w : A — > R 
the correlation function 

p(n) = J A v w o f n dfi — J A v dfi J A w dp, 

decays exponentially: p(n) = 0(r n ) for some r G (0, 1). The corresponding results for 
billiards are complicated by the presence of singularities and unbounded derivatives, 
but these difficulties were eventually overcome by Young (33] and Chernov [5] who 
proved exponential decay of correlations for planar periodic dispersing billiards by 
constructing a certain type of infinite Markov extension called a Young tower [331 El]- 

The situation for flows is considerably more subtle. Whereas mixing Axiom A 
diffeomorphisms always have exponential decay of correlations for Holder observables, 
there exist Axiom A flows that are mixing but at an arbitrarily slow rate [29] [26] . 
Currently, it is not known if exponential decay of correlations is typical for Axiom A 
flows. However, Dolgopyat [2] (see also [12]) proved that typical Axiom A flows are 
rapid mixing (faster than any polynomial rate) for sufficiently regular observables. 
Melbourne [TBI [T9] extended this result to flows with a Poincare map modelled by 
a Young tower, proving for example that the planar periodic Lorentz gas with finite 
horizon is rapid mixing. For this particular flow, Chernov [6] was able to prove 
stretched exponential decay of correlations. 

Thus, generally speaking our understanding of planar periodic dispersing billiards 
and Lorentz gases is now as advanced as in the Axiom A case. However, there are 
a number of classes of planar billiards for which open questions remain [7J [9j (TTJ . In 
this paper, we focus mainly on questions related to decay of correlations for Lorentz 
gas flows, building on [18] (rapid mixing) and [19] (slow mixing). 

Example 1.1 (Dispersing billiards with cusps) Dispersing billiards with cusps, 
where the boundary curves are all dispersing but the interior angles at corner points 
are zero, turned out to be much more involved technically than the usual Lorentz 
gases. For example ergodicity for the dispersing billiard map with cusps was only 
proved in the mid 1990's by Rehacek [27]. By standard arguments, ergodicity implies 
K-mixing 0, Chapter 6] and the Bernoulli property [H [21] for all the hyperbolic 
billiard examples. 

Concerning rates of mixing, Chernov & Markarian [10J proved, using the method 
of Young towers, that decay of correlations for the billiard map satisfies p(n) = 
0((logn) 2 /n). It was anticipated that the logarithmic factor is an artifact of the 
proof, and this factor was removed by Chernov & Zhang [l3j yielding the decay rate 
p(n) = 0(1/ n). It is believed that this result is optimal for the billiard map. This 
rate of decay of correlations is too weak for strong statistical limit laws (see [TJ). 

Since the collisions inside a cusp occur in quick succession, it is reasonable to 
expect a much faster decay rate in continuous time [HI Section 5.6], possibly even 
exponential decay of correlations. Exponential decay seems to be beyond current tech- 
niques even for simpler problems (such as the finite horizon planar periodic Lorentz 
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gas). Previously, no results on decay of correlations for the flow (not even slow mix- 
ing) were available when there are cusps. In this paper we prove rapid mixing for the 
flow. 

A byproduct of our proof is the almost sure invariance principle (ASIP) for the 
flow. We note that the ASIP implies numerous statistical limit laws such as the central 
limit theorem, the law of the iterated logarithm, and their functional versions, see for 
example [251 121]. 

Example 1.2 (Bunimovich flowers) It was discovered by Bunimovich in the 
1970's that billiard tables with focusing boundary components may also show hy- 
perbolic behaviour. For the first examples of such tables, constructed for example in 
[3J, the boundary components are either dispersing, or focusing arcs of circles, subject 
to some further constraints (see Section [3TT1 below). Given their typical shape, such 
billiards are often called Bunimovich flowers. Later the probably best known example 
of a hyperbolic billiard with focusing boundary components, the Bunimovich stadium, 
was introduced in [I]. Focusing hyperbolic billiards have a rich literature by now - 
it is worth pointing out Wojtkowski's contribution, [52] in particular, which clarified 
much of the mechanism behind hyperbolic behaviour. For further historical accounts, 
and detailed proofs of ergodicity, we refer to the monograph [9] and references therein. 

As for rates of mixing, Chernov & Zhang [12] show that in Bunimovich flowers the 
billiard map has decay of correlations p(n) = 0((logn) 3 /n 2 ). (Again, the logarithmic 
factor appears to be an artifact of the proof and it is expected that 1/n 2 is the optimal 
rate.) It follows that the map satisfies statistical limit laws such as the ASIP. See [20] 
for the scalar case and [21] for Revalued observables. 

By [23], the billiard flow immediately inherits the ASIP (and hence the other 
limit laws). However, previous results were unable to establish estimates for decay 
of correlations of the flow. We recall from [T2] that the only effect slowing down the 
decay rate of the billiard map is sliding along a circular arc, where collisions occur 
in quick succession. Just as in the cusp example, it is reasonable to expect that the 
flow mixes better than the map in Bunimovich flowers. 

In this paper, we prove rapid mixing for the flow. 

Example 1.3 (Bunimovich stadia) Bunimovich [1] established hyperbolicity and 
ergodicity for the stadium billiard bounded by two parallel lines connecting two semi- 
circles. Markarian [17] proved decay of correlations 0((logn) 2 /n) for the billiard map, 
and Chernov & Zhang [13] improved this to obtain the optimal rate 1/n. This is too 
weak for strong statistical limit laws; indeed Balint & Gouezel [T] prove a nonstandard 
limit law (nonstandard domain of attraction of the normal distribution with yj n log n 
normalization) for typical observables. The same limit law holds for the flow [Tl [23]. 
In particular, the ASIP fails for both the billiard map and the flow. 

This time, we do not expect the flow to mix more rapidly than the map, but slow 
mixing does not follow from previous results. In this paper, we prove that the stadium 
flow decays at the same slow rate 1/t as the map. 
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Remark 1.4 Dolgopyat's results [14J on rapid mixing for Axiom A flows require 
the observables to be smooth in the flow direction. This restriction is inherited by 
the generalisation in [T8l fl~9] to nonuniformly hyperbolic flows, and by the current 
paper. In particular, our results on rates of mixing for Lorentz gases do not apply 
to certain physically relevant observables such as position and velocity. At present, 
the only result on decay of correlations for Lorentz gas flows that includes general 
Holder observables is Chernov's result for the planar periodic Lorentz gas with finite 
horizon [6]. 

On the other hand, our proof of the ASIP for billiard flows with cusps in Exam- 
ple 11.11 is valid for all Holder observables. 

The remainder of the paper is organised as follows. In Section [2J we recall back- 
ground material on Young towers [33j El], arid on decay of correlations for flows 
possessing a Poincare map modelled by a Young tower [18j [19] . In Section [3j we 
present our new results for Examples 11.11 11.21 and 11.31 

2 Young towers and flows 

Let / : M -> M be a map, and let h : M -> R + be a roof function. Let (j> t : M h -> M h 
be the corresponding suspension flow. Thus the map f(x) = 4>h( x )( x ) is a Poincare 
map for the flow. Conversely, given / and h we define the suspension M h = {(x, u) G 
M x R : < u < h(x)}/ ~ where (x,h(x)) ~ (fx,0). Then the suspension flow is 
given by cj) t (x, u) — (x, u + t) computed modulo identifications. 

In our examples, M is a finite union of Riemannian manifolds. We will say that 
the map / : M — > M is modelled by a Young tower if it satisfies the axioms introduced 
by Young [33], [34]. In particular, there is a set Y C M with Leb(Y) > that possesses 
an appropriate hyperbolic product structure. Furthermore there exists a measurable 
partition of Y into countably many sets Yj, and r : Y — > N with constant values 
ry t on the Y$, such that the induced map F = f r : Y — > Y is smooth and uniformly 
hyperbolic when restricted to any Y i} with respect to the hyperbolic product structure 
of Y. It is important to note that r is not necessarily the first return time to Y, thus 
the corresponding tower map is a Markov extension of the original map / : M — > M. 
For further details of the construction see the original references. 

For examples from billiards, there is a natural invariant measure \x (the Liouville 
measure) that is equivalent to Lebesgue measure on M with L°° density. In what 
follows we always assume that both the discrete return time r (defined on Y) and 
the roof function h (defined on M) are integrable with respect to Lebesgue measure 
and hence \i. The natural flow- invariant probability measure fi h on M h given by 
fi h = fix lebesgue/ J M h dfi coincides with Lebesgue measure for the billiard examples. 

We aim to compute the rate of decay of correlations and statistical limit laws 
for sufficiently regular observables v,w : M h — > R. For m > 0, rj G (0,1), the 
function space C m,ri (M h ) consists of those v : M h — > K for which ||f|| m ,r; = \\ v \\rj + 
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Il^llr) + ' ' ' + 1 1 d™ v \\v < 00 where d t denotes the derivative in the flow direction and 
HI,, = Moo + sup x ^ ummax{h(x)My)}] \v(x,u) - 

We now summarize results from [T8l [19] used in Section [31 Define 

Pv,w(t) = f M h vwo(f) t dfi h - j Mh v dfi h j Mh w dfi h . 

Theorem 2.1 ( |18j ) Let f : M — ► M be a map modelled by a Young tower and 
h : M —>■ IR + be a Holder roof function. Assume 

• Exponential tails: Leb(y G Y : r(y) > n) = 0(7") for some 7 G (0, 1). 

Then typically the flow <p t '■ M h — > M h zs rapid mixing: there exists rj G (0, 1) and 
for any n > 1, there exists m > 0, C > 1, such that 

\pv,w(j^)\ — ||m,r, ||^ 1 1 m.,77 1 j 

for all v, w G C m ' r '(M /l ) and all t > 0. 

Theorem 2.2 ( [19] ) Let f : M — > M 6e a map modelled by a Young tower and 
h : M —> IR + fre a Holder roof function. Assume 

• Polynomial tails: Leb(y G Y : r(y) > n) = 0(l/n^ +1 ) for some (3 > 0, and 

• 1/h G L°° , or more generally there is an integer p > 1 such that l/h p G L°° 
where h p = h + hof + -- - + ho p^ 1 . 

Then typically the flow <pt '■ M h — > M h is polynomial mixing with rate 1/t 13 : there 
exists 7] G (0, 1), m > 0, C > 1, such that 

\Pv,w{t)\ 5: Cll^llm.TjU^llm.r, t , 

for all v, w G C m ' r '(M /l ) and a// i > 0. 

Remark 2.3 In general, Theorems 12. II and 12.21 hold typically subject to a nondegen- 
eracy condition. However, this nondegeneracy condition is automatic for the billiards 
examples considered here because of the contact-like structure [T9] . 

Remark 2.4 Theorem 12. II implies that in situations where Young [33] obtains expo- 
nential decay of correlations for the billiard map / : M — ► M, rapid decay of corre- 
lations holds for the flow <ftf Similarly, Theorem 12.21 implies that in situations where 
Young [34] obtains polynomial decay of correlations for the billiard map / : M — > M 
with rate 0{l/nP), polynomial decay of correlations with rate 0(l/t^) holds for the 

flow 0f 
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Remark 2.5 The condition that h is Holder can be relaxed in Theorems 12.11 and 12.21 

It suffices that h is uniformly piecewise Holder as follows. The billiard examples are 
examples of smooth dynamical systems with at most countably many singularities. 
Write M — (J Mj where the Mj are the maximal subsets of M on which / is smooth. 
It is enough that sup^ H/iIa^Hjj < oo. 

By construction, the Young tower that models the billiard map "respects" the par- 
tition \jMj in the following sense: as F\y { is smooth, |jYj is necessarily a refinement 
of the restriction of |J Mj onto Y. 

Let us make one more remark. Even though this is a very simple observation, 
this is the Key Observation behind the arguments of Section [3J Note that M can 
be viewed as a Poincare cross-section to the flow. It may be the case that there is 
an alternative cross-section M with Poincare map / : M — > M (and roof function 
h : M — *> M + ) that can also be modelled by a Young tower and such that the tail 
distribution of the Young tower for / : M — > M allows for a better control of the 
statistical properties. Note that the corresponding suspension flow <p t : M h — > M h is 
identical to the original flow. 

Statistical limit laws 

Theorem 2.6 ( |21j ) Let f : M — > M be a map modelled by a Young tower and 
h : M — > 1R + be a Holder roof function. Assume 

• Summable tails: Leb(y G Y : r(y) > n) = 0(1 /n 13 ^) for some (3 > 1. 

Then the (vector-valued) ASIP holds: there exists A > 0, and for any Holder 
observable v : M h — > IR d with mean zero there is a d-dimensional Brownian motion 
W , such that (on a possibly enriched probability space) 

J^ v o (j) t dt = W(T) + 0(T 1 2- x ) asT^ oo a.e. 

Remark 2.7 Again, it suffices that the roof function h is uniformly piecewise Holder 
as in Remark 12.51 Similarly, the Holder regularity requirement for the observable v in 
Theorem (221 can be relaxed. Define V(x) = vo<f> t dt. It suffices that V : M — ► M. d 
is uniformly piecewise Holder. 

Remark 2.8 A more difficult question is the ASIP for the time-one map of a flow. 
A general argument of [151 122] yields the scalar ASIP (for real-valued observables) 
in situations where (i) summable decay of correlations is established by certain tech- 
niques, (ii) the class of dynamical systems is closed under time reversal. This includes 
time-one maps for Axiom A flows [22J. Again the result is restricted to observables 
that are sufficiently smooth in the flow direction. 

It is immediate from these considerations that the scalar ASIP holds for time-one 
maps of the flows in the cusp and flower examples (Examples 11.11 and 11.21) . The 
vector-valued ASIP remains an open question even for time-one maps of Axiom A 
flows. 
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3 Billiard flows 



In this section, we specialise to billiard flows and treat the three examples discussed 
in the introduction. We begin with Example 11.21 (flowers) since this is the simplest. 
Then we treat Examples 11.11 (cusps) and 11.31 (stadia). 

Let Q be a compact domain in M 2 with piecewise C 3 boundary dQ and let us 
consider the corresponding billiard dynamics in Q - the motion of a point particle 
that travels with unit speed and bounces off dQ elastically (the angle of incidence 
is equal to the angle of reflection). The resulting flow <p t is a three-dimensional 
volume preserving flow. There is a natural two-dimensional cross-section M = dQ x 
[— 7r/2,7r/2] corresponding to collisions, and the Poincare map / : M — > M is called 
the billiard map. We choose coordinates x = (r, <p) on M, where r (the position) is 
the arclength parameter along dQ and tp (the outgoing velocity) is the angle made 
by the reflected trajectory and the normal to the boundary. The natural invariant 
measure for the billiard map is is the Liouville measure dfi = cos ip dr d(p. 

By components or arcs we mean the maximal connected C 3 -smooth pieces of the 
boundary dQ. We use the same terminology for the corresponding (two-dimensional) 
components of the phase space M. We assume that the curvature of any component 
has a fixed sign, thus we may consider convex inwards (dispersing), convex outwards 
(focusing) and neutral arcs, and decompose the phase space as M = M + U M~ U M°, 
respectively. 

The dynamics of the billiard map/flow may depend sensitively on further geo- 
metrical properties of dQ, including the types of the components present and their 
distances and angles. Below we restrict to certain physically relevant examples and 
cite the literature for their known dynamical properties. For further details on billiard 
dynamics in general, see [9J and references therein. 

Our analysis relies on the standard fact that the billiard map / and the roof 
function (collision time) h are piecewise Holder continuous. A simple argument shows 
that the Holder exponent is \ for both / and h. Note that h(x) = \f(x) — x\ so it 
suffices to verify that / is ^-Holder. 

It is clear that / is smooth except near preimages of tangencies (grazing collisions) 
and writing / = (f r , ftp) a standard calculation 0, Section 2.11] shows that /' is 
bounded except for singularities where it behaves like 1/ cos(f ip (x)). Writing ip = 
n/2 — (p and suppressing the r-coordinates, we have fL(x) ~ 1/ sm(f^(x)) ~ 1/ (f^(x)) 
so that ((fip(x)) 2 )' ~ 1. Hence fy(x) ~ ip 1 ^ 2 which is ^-Holder as required. 

3.1 Bunimovich flowers 

Bunimovich 0] studied hyperbolicity and ergodicity for a class of billiard tables Q C 
M 2 such that (i) M° = (no neutral components); the dispersing components M + 
may have arbitrary geometry, in contrast (ii) M~ consists of circular arcs all strictly 
smaller than a semicircle; (iii) if such a circular arc was extended to a full circle, this 
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circle would be contained in Q; (iv) the neighbouring components are transversal (no 
cusps). Actually, the conditions can be relaxed and even neutral components can be 
allowed as long as certain technical assumptions are satisfied, see [12] for details. 

Chernov & Zhang [12] show that for Bunimovich flowers the billiard map / : 
M — > M can be modelled by a Young tower with tail estimate roughly 0(l/n 3 ). It 
then follows from Young [M] that the map has decay of correlations roughly 0(1/ n 2 ). 
By [2TJ, the vector- valued ASIP holds for Holder observables for both the billiard 
map and the billiard flow. 

However, there is no immediate conclusion for the decay of correlations of the flow. 
(Exponential tails are required in Theorem 12.11 while Theorem 12. 21 would yield at best 
0(l/t 2 ) decay. In any case, the assumption that h is bounded below in Theorem 12.21 
is violated.) Nevertheless, we now prove rapid mixing for the flow for sufficiently 
smooth observables. (The scalar ASIP for the time-one map for the flow then follows 
from Remark 12. 81 ) 

Theorem 3.1 Let <p t be the flow corresponding to a Bunimovich flower. Then <p t is 
rapid mixing (in the sense of Theorem \2.1\) . 

Proof The argument in Chernov & Zhang [121 p. 1546] demonstrates the existence 
of an alternative cross-section M C M (with Poincare map / : M — > M) that is 
modelled by a Young^ tower with exponential tails. We claim that the corresponding 
collision time h : M — > M + satisfies the uniformly piecewise Holder condition in 
Remark 12.51 expressed of course in terms of M, f, h. Rapid mixing follows from 
Theorem 12.11 

First, we recall the definition of M in [12]. As mentioned above we write M = 
M + U M~ where M + corresponds to dispersing arcs and M~ corresponds to focusing 
arcs. Then M is given by M = M + U E where E C M~ consists of only the first 
(sliding) collisions at each focusing arc, so E = M~ fl f(M + ). 

Next we verify that h G L°°. Since h = h on dispersing arcs we can restrict 
attention to a single focusing arc T(c M~), and in particular the single first collision 
set Er = T H E. Consider x G Er- There is an integer n > 2 such that f^x G T for 
j = 1, . . . ,n and f n+l x G" T. Hence 

h = g + hof n 1 (3.1) 

where g(x) = h(x) + h(fx) + • • • + h(f n ~ 1 x) is the amount of time it takes to "slide" 
along the arc T. Hence |/ils|c» < |T| + |/i|oo and so h is bounded. 

Finally, we verify that there is a uniform Holder constant for h on each partition 
element. Again we may restrict to a single first collision set Er C T where T is a 
focusing arc, but this time it is necessary to consider the finer partition Er = Un>i 
where E n = {x G E : f^x G T for j = 1, . . . ,n and/ n+1 x G^ T}. (In other words, 
/ = f n+1 on E n .) Note that the partition elements E n coincide with the smoothness 
components of / : M — > M, which is in agreement with Remark 12.51 
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Let R denote the radius of the focusing arc Y. Let (r, tp) denote the standard 
coordinates and note that sliding occurs for ip ~ tt/2 so it is convenient to introduce 
ip — 7r/2 — (p. It is immediate from the geometry of the circle that 

f(x) = (r + 2Rip,ip), h(x) = R sin i/j 

for all x — (r, ip) G Y n / -1 r. Hence 

|<7(#i) — g{x 2 )\ < sin?/'! — sin ^2 1 < nir^i — ^2 1 (3.2) 

for all Xi, x 2 G E n . On the other hand, since the angle ip remains constant during the 
n grazing iterations, it is elementary that there is a constant C\ depending only on 
r such that Ci/(n + 1) < ip < d/n for all x G E n . Hence |^i -ih\<C x (±- ^) < 
Ci/n 2 and so 

hh-^l'^Cfn- 1 . (3.3) 

Combining estimates (I3.2p and (13. 3p . we obtain \g{x\) — g(x 2 )\ < C 2 \ipi — ip 2 \^ where 
1 

C 2 = RC{ depends only on the arc Y . 

In addition, since f n (r, ■0) = (r + 2nRip,ip) on i? n and ft, ^-Holder, \h{f n x\) — 

Hf n x 2 )\ <C 3 (\r 1 -r 2 \ i +(2Rn + l)\t(; 1 -H^ < ^(h - r 2 \ + |Vi - ^1*)* < 
C^(\ri — r 2 \ 2 + — ^2 1*)- Combining the estimates for g and hof n , we have shown 
that /i is Holder on with constant independent of n (and exponent |). Since there 
are only finitely many arcs, this completes the proof. ■ 

3.2 Dispersing billiards with cusps 

Chernov & Markarian [10] studied dispersing billiards with cusps, where the boundary 
curves are all dispersing - that is, M = M + - but the interior angles at corner points 
are zero. By [TQl EE], the billiard map / : M — > M can be modelled by a Young 
tower with tail estimate 0(l/n 2 ). It follows from Young [3lj that the map has decay 
of correlations 0(l/n). This is too weak for strong statistical limit laws (see pQ). 
Nevertheless we now prove rapid mixing and the ASIP for the flow. 

Theorem 3.2 Let <j) t be the flow corresponding to a billiard table with cusps. Then 
<fit is rapid mixing (in the sense of Theorem \2.1\) and satisfies the ASIP (in the sense 
of Theorem \2.6\) . 

Proof Following Chernov & Markarian [10] we define M by excluding a neighbour- 
hood of each cusp. The new collision map / : M — > M is modelled by a Young tower 
with exponential tails |10j . By Theorems 12.11 and 12. 6[ it again suffices to show that 
the new collision time h : M — > M + satisfies the uniformly piecewise Holder condition 
in Remark 12.51 
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Consider a single cusp and let E n be the set of those points in M that spend 
exactly n iterates in the cusp before returning to M. Note again that, in accordance 
with Remark 12. 5[ the sets E n coincide with the smoothness components of /. The 
calculation in [T0~1 p. 738] shows that h is bounded on E n independent of n. Hence h 
is bounded on the whole of M. 

It remains to find a uniform Holder constant for h on each E n . Explicit calculations 
are more difficult than in Section 13.11 so we search for coarser estimates. We claim 
that there are positive constants «i, a 2 , « 3 , C*i, C 2 (independent of n) such that 

(i) \h(xi) — h(x 2 )\ < Cin ai \xi — x 2 \ a2 , and 

(ii) \h( Xl ) -h{x 2 )\ < C 2 n- a \ 

for all xi, x 2 G E n . It then follows that \h(xi) — h(x 2 )\ 1+1 < C\x\ — x 2 \ a21 with 
7 = atz/a.i, C = CjC 2 , and so h is uniformly piecewise Holder as required, with 
exponent 027/(7 + !)• 

We verify the claim with a\ = ^r, a 2 = \ and a 3 = | (and hence Holder exponent 
i). Throughout, C is a uniform constant that can change from line to line. Let 
Xi,x 2 G E n . Recall that / and h are Holder with exponent |. Furthermore, by [TU1 
Proposition 4.1], 

\Pxi ~ f j x 2 \ < Cn 5 / 3 \f Xl - fx 2 \ < Cn 5 ^\ Xl - x 2 \ 1/2 

for j — 1, . . . , n. We deduce that 

- h{x 2 )\ < \h{ Xl ) - h(x 2 )\ + Hfxx) - h{Px 2 )\ 

3=1 
n 

<C\x 1 -x 2 \ 1 / 2 + Cj2\rx 1 -fx 2 \ 1 / 2 

< C\x x - x 2 \ l/2 + nCn^\x x - x 2 | 1/4 < Cn 11/6 |n - x 2 \ 1/4 

establishing estimate (i). 

By pH P- 738, last line], we have h(f Xj ) + ■■■ + h{f n - l Xj) < Cn' 1 for j = 1, 2. 
As shown in [TOf pp. 748-749], the cell E n has diameter of order n~ 2 ^ 3 so that 
|/i(a?x) — h(x 2 )\ < C\x\ — x^ 1 / 2 < CrT 1 ! 3 . By time-reversibility of the construc- 
tion, \h(f n xi) — h(f n x 2 )\ < Cn^ 1 / 3 . Altogether, we have \h(xi) — h(x 2 )\ < Cn^ 1 ^ 3 
establishing estimate (ii). (The estimates of [10] that we have used are established in 
the first instance for a special billiard table with three cusps, but then extended to 
the general situation in [TUl Section 6].) I 

Remark 3.3 Again, we stress that the vector-valued ASIP for the flow holds for gen- 
eral (piecewise) Holder observables, whereas rapid mixing is restricted to sufficiently 
smooth observables (as is the scalar ASIP for the time-one map of the flow which 
holds by Remark 12.81) . 
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3.3 Bunimovich stadia 



Bunimovich [4] established hyperbolicity and ergodicity for the stadium billiard 
bounded by two parallel lines (M°) connecting two semicircles (M~). By [T71 113] , 
the billiard map / : M — > M can be modelled by a Young tower with tail estimate 
0(l/n 2 ). It follows from Young [34] that the map has decay of correlations 0(l/n) for 
Holder observables. This is too weak for strong statistical limit laws; indeed Balint & 
Gouezel pQ prove a nonstandard limit law (nonstandard domain of attraction of the 
normal distribution with \J n log n normalization) for typical observables. By |23j, the 
same limit law holds for the flow. In particular, the ASIP fails for both the billiard 
map and the flow. 

This time, we do not expect the flow to mix more rapidly than the map. Nor 
can we apply Theorem 12.21 directly since the roof function is not bounded below. 
Nevertheless the conclusion of Theorem 12.21 is valid. 

Theorem 3.4 Let <p t be the flow corresponding to a Bunimovich stadium. Then (f) t 
is polynomially mixing with rate 1/t. 

Proof The argument in Chernov & Zhang [12[ p. 1548] demonstrates the existence 
of an alternative cross-section M C M (with Poincare map / : M — > M) that is 
modelled by a Young tower with exponential tails. However, the corresponding roof 
function h : M — > R + is unbounded, so Theorem 12.11 does not apply this time. As 
mentioned above, they show that the corresponding tower for / : M — *■ M has tails 
decaying at the rate 0(l/n 2 ). We note that the tower for M is built over the same 
base Y as the tower for M but is strictly higher (since it incorporates the returns 
from M to M). 

Whereas the cross-section M in [12] takes account of both sliding along circular 
arcs in M~ and bouncing between the parallel straight lines in M°, we define an 
intermediate cross-section M' with M C M' C M that takes only account of sliding. 
That is, we define M' = M° U E where E — M~ n/(Af°) consists of the first (sliding) 
collisions in M~ . 

It is immediate that the tower for /' : M' M' has tail decay rate no worse than 
the tower for it shares the same base Y as the other two towers, but is lower than 
the tower for M (and higher than the tower for M). In particular the tower for M' 
has tail decay rate 0(l/n 2 ). 

Let h! : M' — ► R + denote the new roof function. It remains to verify the hypothe- 
ses on h! in Theorem 12.21 The proof that h! is uniformly piecewise Holder is identical 
to that in Section 13.11 We claim that h' is bounded below for p = 2. Indeed, h! 
can approach zero only at points x G M° which are close to one of the endpoints 
of a straight boundary component, that is, to the corner made with M~, and have 
velocity almost tangent to the straight boundary component. But then fx is about 
to undergo a long sequence of sliding collisions so that h'(fx) has magnitude approx- 
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imately the length of a semi-circular arc in M . Hence h' 2 is bounded away from zero 
as required. I 
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